Origin of the Mott Gap 
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We show exactly that the only charged excitations that exist in the strong-coupling limit of the 
half-filled Hubbard model are gapped composite excitations generated by the dynamics of the charge 
2e boson that appears upon explicit integration of the high-energy scale. At every momentum, such 
excitations have non-zero spectral weight at two distinct energy scales separated by the on-site 
repulsion U. The result is a gap in the spectrum for the composite excitations accompanied by 
a discontinuous vanishing of the density of states at the chemical potential when U exceeds the 
bandwidth. Consequently, we resolve the long-standing problem of the cause of the charge gap in a 
half-filled band in the absence of symmetry breaking. 



I. INTRODUCTION 

In 1949, Sir Neville Motti proposed that transition 
metal oxides with half-filled bands possess a gap in the 
single-particle spectrum that is due entirely to the en- 
ergy cost for placing two electrons on the same site. This 
explanation is clearly incomplete because even in the sim- 
plest model of a Mott insulator, the Hubbard model, none 
of the eigenstates have definite local occupation. Conse- 
quently, the charge gap in transition metal oxides does 
not have the simple interpretation as the energy gap be- 
tween bands that represent electron motion among singly 
and doubly-occupied sites. What then are the degrees of 
freedom that are being gapped? Because mobile doubly 
occupied sites would be inconsistent with an insulating 
state, some^ have argued that in a Mott insulator, dou- 
ble occupancy is localized whereas in the metal doubly 
occupied sites form an itinerant fluid. Such localization 
requires a dynamical degree of freedom which has not 
been ennunciated despite numerous simulations which 
display a Mott gap^. In fact, the origin of the dynamical 
degree of freedom that generates the elementary excita- 
tions responsible for the Mott gap is the essential problem 
of Mottness. Knowledge of this degree of freedom and 
the excitations it mediates are crucial to the physics of 
high-temperature copper-oxide superconductors as they 
are doped Mott insulators. Indeed, the extreme difficulty 
in unearthing the mechanism for the localization of dou- 
ble occupancy led Laughlin 5 to suggest that charge gaps 
in homogeneous time-resersal systems are impossible. 

In this paper, we construct explicitly the dynamical 
degrees of freedom that account for the Mott gap in the 
absence of any symmetry breaking. There are two key 
elements to our proof. First, we show that the exact low- 
energy theory for a half-filled band described by the Hub- 
bard model has no bare propagating degrees of freedom 
of any kind. Second, we identify two dispersing degrees of 
freedom or composite excitations which lead to a turn-on 
of the spectral weight centered at ±U/2. If U > 8t, the 
composite excitations are orthogonal to one another in 
that the spectral weight they produce never exists in the 
same energy range. The result is a gap in the spectrum. 
In terms of the UV variables, the composite excitations 



represent bound states^ involving double occupancy or 
double holes and are the fundamental excitations that 
define the lower and upper Hubbard bands. 



II. LOW-ENERGY ACTION 

As pointed out by Laughlin£, no one has identified the 
band structure of the elementary particles whose spec- 
trum becomes gapped at half-filling. In this paper, we 
show that this can be done by utilizing the methods^ we 
have recently developed to explicitly integrate out the 
degrees of freedom far from the chemical potential in the 
Hubbard model. We consider the Hubbard model on a 
square lattice in the limit in which the bands are well 
separated, that is, the on-site interaction U exceeds the 
bandwidth, W — St, t the hopping matrix element. At 
half-filling, the chemical potential lies in the Mott gap. 
As a consequence, both the degrees of freedom above 
and below the chemical potential must be integrated out 
if one wishes to construct a low-energy theory of the 
Mott insulator. This can be done by introducing^ two 
new fermionic fields which when constrained appropri- 
ately will correspond to the creation of double occupancy, 
Di, and double holes, Di. In Lorentzian signature, the 
Lagrangian which makes this integration possible, 
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iD ] t) - iDD - —{D^D - bb^) 



-0bD + h.c. + s6cpl(D - 0c ]Ci ) 



+ s6(p\b - 9c\c\) + h.c. 



(1) 



contains the two constraint charge ±2e bosonic fields, (p\ 
(charge 2e) and (p\ (charge — 2e) which enter the the- 
ory as Lagrange multipliers for the creation of double 
occupancy and double holes, respectively. Mathemati- 
cally, they are analogous to a in the non-linear sigma 
model. All operators in Eq. (pj have the same site in- 
dex which is summed over. The Lagrangian also contains 
the integration, d 2 0, over the complex Grassman, 6, and 
bi = Ylj 9ij c i.<jVa c j,-a is a bond-singlet operator where 
ct creates a fermion on site i with spin a, gij = 1 iff i 
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and j are nearest neighbours, V[ = — Vi = 1 and s and 
s are constants appearing in the constraint which have 
units of energy. The theory (fTJ) is completely equivalent 
to the Hubbard model. That this is so can be seen by 
integrating out ipi and <fi followed immediately by an 
integation over Di and Di in the partition function, 

Z = J[Dc 2?c t VD VD*< Vlp exp~ fo L " vdt . (2) 

The ip and <p integrations (over the real and imaginary 
parts) are precisely a representation of (a series of) 5- 
f unctions of the form, 



J d9D t - J d6 6c iA Ci^ , 



and 



dOD; - 



d8 04Al 



(3) 



(4) 



We must now integrate over the Di and and Di. The 
dynamical terms yield, 



1 C i>jCi<7- 



(5) 



The terms proportional to U lead to 

U 
~2 



d 2 e 



= U 



(6) 



the standard interaction term in the Hubbard model in 
the Lorentzian signature. Finally, the terms proportional 
to V a yield 



d 2 e 



•E 



h.c, 



h.c. 



(J) 



after the if and Di integrations, whereas the (p and Di 
integrations yield the same final result except rij^ a is re- 
placed by (1— nj— a ). Hence, these terms add together to 
generate the kinetic term in the Hubbard model. Adding 
all these results together leads to / d 2 998Cnubb = Aiubb, 
the Lagrangian for the Hubbard model. Hence, Eq. (fTJ) 
is an equivalent way of writing the Hubbard model in 
which two canonical fermions describe the physics on the 
[/—scale. 

Since the physics on the [/—scale has been cleanly iden- 
tified, the exact low-energy theory, 
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(s</? f - s*ip) c T c x + h.c, 
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can be constructed by explicitly integrating out the mas- 



sive fields Di and D^. Here, L± 



di 



- 1 - 2 ' 



This 



tegration is straightfoward as it is strictly Gaussian. It 
is important to appreciate that the resulting theory is 
exact, and completely equivalent to the Hubbard model. 

We would like to understand the physics of this model. 
It is clear from this presentation that the bosonic and 
fcrmionic fields do not represent weakly coupled degrees 
of freedom, as there are no quadratic terms in and the 
would-be boson propagators have no poles. Presumably, 
one would like to find an appropriate continuum limit, 
but the identification of the correct continuum limit is 
difficult. What we will find here is that Eq. ([5]) is a much 
better theory to work with than the Hubbard model, 
as it contains the seeds of the degrees of freedom that 
emerge at low energies. The Hubbard model contains 
only strongly interacting electrons, and any continuum 
limit that might be considered would presumably miss 
any collective degrees of freedom. 

To proceed, we begin by recalling what happens in the 
free fermion Landau theory. There, the continuum limit 
is trivial to take, as we just scale towards the free fermion 
UV fixed point. The appropriate renormalization group 
flow is obtained^ by scaling momenta towards the Fermi 
surface, k = kp + 1, 1 — ► 0. The fundamental reason that 
the latter is done is that that is where the spectral density 
lies. The correct degrees of freedom in the IR give rise 
to this spectral density. At the level of the Lagrangian, 
the spectral density is determined by the vanishing of the 
coefficient of the quadratic terms. The spectral density is 
highly peaked and the effects of renormalization are only 
to give weakly interacting (dressed) fermions. 

We emulate this approach by determining where Eq. 
(0) predicts the spectral density to lie. Fortunately, we 
will find that this is highly peaked (at what we will call 
the upper and lower Hubbard bands), and so one might 
hope that the scaling towards that locus is well-defined. 
We will not go so far in this paper as to claim that a 
continuum limit exists, but we will use this insight to 
establish that a dynamical Mott gap emerges. The root 
cause of this effect is that hidden in Eq. (JSJ is the dynam- 
ical degree of freedom that mediates the Mott gap: both 
the electrons and bosons are locked into bound states 
and hence cannot propagate independently. As a result 
the most relevant term in the Lagrangian arises from 
the boson-fermion interaction. This is a purely strongly 
coupled effect which makes Mottness analogous to other 
problems of strong interactions, for example QCD. An 
important difference with QCD which Eq. |(5J) lays plain 
is that for Mottness, the exact low-energy Lagrangian 
may be derived. This should enable an identification of 
the proper collective degrees of freedom. 

To proceed, we switch to frequency and momentum 
space and specialize to a square lattice as our focus 
is the copper-oxide plane of the cuprates. Defining 



ip(t) = J duj 



(full the energy dispersion, e 



(k) 



4 J^fi cos{k^a/2) cos(p M a), where k and p are the center 
(8) of mass and relative momenta of the fermion pair, and 
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the Fourier transform of bi, 



we arrive at the exact working expression, 



(9) 



rhf _ 
^IR — 



S| t \ s \ ~t 



(w - {7/2) 



(w + {7/2) 



{7 2 -4w 



x (ck/2+ P ,TC k /2_ Pii ) w + h.c. 



(10) 



for the low-energy Lagrangian where we have suppressed 
the implied integration over frequency and introduced 
the coupling constants, 



(k) -U + te^ + 2uj 



« P k) H = 



U-2w 
U + te ( p } + 2lu 
U + 2uj 



(11) 



which play a central role in this theory. They, in fact, 
will determine the spectral weight in the lower Hubbard 
(LHB) and upper Hubbard (UHB) bands, respectively. 
Note that in all of these expressions, to is the frequency 
of the boson field ip or (p. As we have retained the full 
frequency dependence, we will be able to determine the 
complete dynamics. 

If the bosons were weakly coupled propagating degrees 
of freedom, setting the coefficient of the quadratic terms 
(in the Lagrangian) to zero would determine their disper- 
sion. However, the coefficients of the naively quadratic 
terms never vanish for any momentum and frequency. 
Hence, on the surface of it, neither the bosons nor the 
electrons propagate and the spectral weight vanishes at 
all energies. The correct theory of the Mott gap should 
yield, however, a non-zero spectral weight in the UHB 
and LHB. Identifying this physics requires that we re- 
evaluate what should properly be considered to be a ki- 
netic term. The structure of the frequency and momen- 
tum dependence of Eq. (fT0|) suggests that the opera- 
tors ip^cc and (pec play a central role and they deter- 
mine where the spectral weight resides. These operators 
might then be thought of as the kinetic terms for com- 
posite excitations mediated by the charge ±2e bosonic 
fields (loosely speaking, we might think of this as occur- 
ring because of the formation of bound states) . Such an 
interpretation is warranted because the spin-spin interac- 
tion and all higher-order operators contained in the \b\ 2 
term are at least proportional to a 4 and hence are all 
sub-dominant to the composite interaction terms. Con- 
sequently, at the level of the Lagrangian, the turn-on of 
the spectral weight is governed by the vanishing of the 
coefficients of the coupled boson-fermion terms. 

That novel dynamics emerges from Eq. (JTUJ) can be 
seen by inspection of the coefficients (TTT|) . We note 
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FIG. 1: a) Diamond-shaped surface in momentum space 
where the particle dispersion changes sign, b) Turn-on of the 
spectral weight in the upper and lower Hubbard bands for the 
composite excitations as a function of energy and momentum. 
In the UHB, the spectral density is determined to 7 P while 
for the LHB it is governed by y p . The corresponding oper- 
ators which describe the turn-on of the spectral weight are 
the composite excitations ^ cc (UHB) and (pec (LHB). The 
electron spectral density is determined by an overlap (see Eq. 
(|16[) ) with these propagating collective modes. 



that the frequency poles appearing in the various terms 
of the Lagrangian are an artifact of our normaliza- 
tion, and could be absorbed into a redefinition of fields: 
ip u — > — 2uj/U ip u , <p u — ► y/l + 2uj/U <p u , and 
(cc) w — ► y/l — Alu 2 /U 2 (cc) w . These scalings recast the 
Lagrangian as 



L 



hf 
IR 



2 >. 



1 U 



<P- 



t 2 



(12) 



+S7 p k) (w)^ jk C k /2 + p^/2+ w ',TCk/2-p, W /2-a ) M 
+ S*7p k) ( w )^-w,kCk/2+p,w/2+ W ',T c k/2-p,( i ;/2-t J '4 

(13) 



h.c. 



Effectively, we have rescaled the coefficients a, a to 

rww \ -U + te i v ) +2uj , — 

7 P k) M = y/l + 2w/U 

U + te p k) + 2lu 



7 P k) M = r 



VI - 2u/U, (14) 



while the coefficients of other terms in the Lagrangian 
are just constants. The first thing to notice about these 
expressions is that the boson frequency appears in the 
combinations U =F 2u>. What this will ultimately mean is 
that the analytic structure is concentrated around lo = 
±{7/2. 



A. Propagating Degrees of Freedom at Half-filling: 
Mott Gap 

To determine where the spectral weight resides, we cal- 
culate where the coefficients 7p and 7p vanish. Consider 

initially k = so that the dispersion simplifies to £ p = 
4 J2 U cos a Pv When ui = ±U /2, 7 p (jp) vanish along the 
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momentum surface denned by X^ C0S fV = 0- This sur- 
face corresponds to the diamond ap = (ap x , ±7r — ap x ) 
depicted in Fig. (QJi). These features define the center of 
the LHB (-U/2) and UHB (U/2) for the composite exci- 
tations. For all momenta outside the diamond, |p| > tt, 
7 p vanishes for U/2 < uj < U/2 + At while 7 p = for 
-U/2 < uj < -U/2 + it. Within the diamond, |p| < tt, 
7 P = in the energy range U/2 — At < uj < U/2 while in 
the interval [— U/2 — At, —U /2] the coefficient 7 P vanishes. 
Consequently, for each momentum, spectral weight turns 
on at two distinct energies, one in the LHB (7 P = 0) and 
the other in the UHB (7 P = 0) with a separation of U . 
Note that 7 P and 7 P never vanish at the same energy pro- 
vided that U > 8t. Consequently, for U > W, a hard gap 
opens (the Mott gap) in the spectrum and the excitations 
defined by the vanishing of 7 and 7 propagate indepen- 
dently above and below the gap, respectively. In terms 
of the composite excitations, the Mott gap opens contin- 
uously but the spectral weight at the chemical potential 
rises discontinuously as is seen in numerical calculations 
on finite-dimensional lattice s 11 ' 12 but in contrast to the 
d = oc 4 solution. The composite excitations which lead 
to the turn-on of the spectral weight correspond to the 
bound states of ip^cc (UHB) and (pec (LHB); they rep- 
resent the collective modes or in essence the propagat- 
ing charge degrees of freedom of the half-filled Hubbard 
model. This is our principal conclusion. As our analysis 
thus far is exact, we conclude that in the absence of any 
symmetry breaking, the coefficients 7 P and 7 P determine 
the dispersion for the excitations that comprise the here- 
to-fore undefined 5 UHB and LHB. The center-of-mass 
momentum k simply shifts the momentum at which 
changes sign, thereby keeping the Mott gap intact. 

Ultimately, it is the overlap between the composite ex- 
citations and the bare electrons that determines the turn- 
on of the electron spectral density. Consequently, the gap 
in the electron spectrum is at least that of the composite 
excitations. To determine the overlap, it is tempting to 
complete the square on the ip'cc term bringing it into a 
quadratic form, with 'J = Aip + Bcc. This would 

lead to composite excitations having charge 2e, a vanish- 
ing of the overlap and hence no electron spectral density 
of any kind. However, the actual excitations that under- 
lie the operator ip^cc correspond to a linear combination 
of charge e objects, c' and (p'c. In terms of the UV vari- 
ables, the latter can be thought of as a doubly occupied 
site bound to a hole. To support this claim, we construct 
the exact representation of the electron creation operator 
at low energies. This can be done by adding to the start- 
ing Lagrangian a source term that couples to the current, 
Jj, that generates the canonical electron operator when 
the constraint is solved. In this case, 



creation operator in the IR at half-filling 



rhf 



iuv + / d 2 9J t , a V a D\ Cl ^ a e + VaOci^Di + h.c. 



is the correct transformation to generate the canonical 
electron operator in the UV. If we now integrate the par- 
tition function over Di and D;, we find that the electron 



S( Pi 



bio. 



(15) 



is indeed a sum of two composite excitations, the first 
having to do with spin fluctuations (b'c) and the other 
with high-energy physics, ip^c and ipe, that is, excita- 
tions in the UHB and LHB, respectively. It is important 
to note that Eq. (fT5"|) is the exact expression for the 
low-energy electron at half-filling. Consequently, we for- 
mulate the overlap 



= |(c t |c t )(c t |* t )| 2 P* 



(16) 



for the the physical process of passing an electron through 
a Mott insulator in terms of the overlap between the bare 
electron with the low-energy excitations of Eq. (|15p . 
(c|c), and the overlap with the propagating degrees of 
freedom, (c\^) with Py, the propagator for the compos- 
ite excitations. As a result of the dependence on the 
bosonic fields in Eq. (fT5")) . O contains desructive inter- 
ference between states above and below the chemical po- 
tential. Such destructive interference between excitations 
across the chemical potential leads to a vanishing of the 
spectral weight at low energie s 11 ' 13 ' 14 . Consequently, the 
turn-on of the electron spectral weight cannot be viewed 
simply as a sum of the spectral weight for the composite 
excitations. As a result of the destructive interference, 
the gap in the electron spectrum will always exceed that 
for the composite excitations. Hence, establishing (Fig. 
(H|)) that the composite excitations display a gap is a suf- 
ficient condition for the existence of a charge gap in the 
electron spectrum, a key conclusion of this work. 



B. Electron Spectral Function 

We confirm the argument in the previous section on the 
origin of the gap in the electron basis by an explicit calcu- 
lation of the electron spectral function. Because the ac- 
tion lacks any derivative terms with respect to ipi, we can 
treat ip as a spatially homogeneous field. While A priori, 
such gradient terms with respect to ipi are possible, their 
presence at half-filling would indicate that freely propa- 
gating bosonic degree of freedom exist at half-filling. The 
absence of such terms at half-filling makes it possible to 
identify that the only propagating degrees of freedom at 
half-filling are gapped composite excitations. 

We proceed by rewriting the coefficient of the boson- 
fcrmi terms as 



A(k,u,(p,4>) = -s(«^ t - (p) 



+ ( 



st 



st 



U-2uj-i5 Lf ' U + 2uj + i5 



and 



a(k) = 2t(cos(k x ) + cos(fc y )). 



(p)a{k) 



(17) 
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Any non-trivial dynamics underlying the Mott gap will 
arise only from the the second term in A(k,u>,ip,<p). 
Upon Wick rotation, ip — > ip and ip* — > ip*, we rewrite 
the single-particle electron Green function as 

G(k,uj)= { dip { d<pG(k,uj,ip,<p>)exp~ 1 ' (18) 



where Cyiott is the IR Lagrangian with the \b\ 2 term 
dropped and 



\A(k,u;,p,p)\ 2 +i5 



(19) 



As our analysis thus far demonstrates that the \b\ 2 term 
has no bearing on the Mott gap justifies our use of the 
truncated Lagrangian, £Mott which has only the charge 
degrees of freedom, because of the i 6 in the gap function, 
A(k,u>, p, 0), the imaginary part of the Green function 

%G(k,uj,p,p) = lim [(U - 2u) 2 + S 2 ] \(U + 2uj) 2 + S 2 ] 



A 2 + {2A(<p + >p) + B 2 ) S 2 + 0(6 4 ) 
(U -2uo) 2 {U + 2lu) 2 



B 



-6(A) 



(20) 



is explicitly non-zero. We have defined 

A = [U 2 - 4lu 2 - 2a k (U + 2uj)] ip 
+ [U 2 - Au: 2 - 2a k (U - 2uj)] p 

2ip(2u + a k ) + 2p{2uj + a k ). (21) 



B 



Note the arguments of the 5 functions are closely related 
to the coefficients 7 and 7 that led to the turn-on of the 
spectral weight. However, in the electron basis, the spec- 
tral weight is not a simple sum of the spectral weight at 
±{7/2. As a result of the integration over ifi and (pi, the 
spectral function for the electrons arises from a compli- 
cated interference between excitations in the LHB and 
UHB. Consequently, to complete the calculation, we per- 
formed the (pi and <p integrations numerically. The resul- 
tant electron spectral function for U = 8t shown in Fig. 
@ demonstrates clearly that a Mott gap exists and the 
spectral weight is momentum dependent. This calcula- 
tion supports the physical argument made in the previ- 
ous section that the Mott gap in the electron basis arises 
from a non-trivial interference between the excitations 
at ±{7/2. At (7r,7r), the spectral weight lies predomi- 
nantly in the UHB whereas at (0, 0) it lies in the LHB. 
Consequently, the real part of the Green function must 
change sign** along some momentum surface that lies be- 
tween these these two extreme momenta. The location 
of the zero surface or Luttinger surface is the Fermi sur- 
face of the non-interacting system as it must b o 16 ' 17 for 
the half-filled system with particle-hole symmetry. We 
find then that the Mott gap arises from the dynamics of 
the two charge \2e\ bosonic fields. This is the first time 
the Mott gap has been derived dynamically, in partic- 
ular by a collective degree of freedom of the lower and 
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FIG. 2: Electron spectral function corresponding to Lmou for 
U = St. The gap here is generated entirely from the dynamics 
of the charge 2e bosonic fields that emerge from integrating 
out the upper and lower Hubbard bands at half-filling. 



upper Hubbard bands. Relative to the gap in the spec- 
trum for the composite excitations that diagonalise the 
fcrmion-boson terms in Eq. (fT5|) , the gap in the electron 
spectrum is larger. This is not surprising as the bare 
electrons do not have unit overlap with the composite 
excitations. While our treatment of the charge ±2e bo- 
son is approximate, it does suffice to capture the essence 
of of the collective mode, namely it mixes all sectors with 
varying numbers of doubly occupied sites. In addition, 
we anticipate that the electron spectral function should 
evolve as the Mott transition is approached in a similar 
fashion to that in terms of the composite particle basis. 
That is, the gap should close continuously without a co- 
herence peak at zero energy. Conseqently, the spectral 
weight at the chemical potential should jump discontinu- 
ously from zero to the value in the free system at the Mott 
transition as is seen in simulations of the Mott transition 
in finite-dimensional systems. In addition, the momen- 
tum dependence of the spectral function is identical to 
that obtained by dynamical mean-field calculation ;] 18 ! 19 
thereby lending crecedence to such cluster calculations^ 
near the Mott transition. 



III. FINAL REMARKS 

What this analysis demonstrates is that the spin-spin 
interaction, contained in the \b\ 2 term, plays a spectator 
role in the generation of the Mott gap. Nonetheless, there 
is a natural candidate for the antiferromagnetic order, 
namely Bij = {9ijf\ci,^c-j,i) ■ The vacuum expectation 
value of this quantity is clearly non-zero as it is easily ob- 
tained from a functional derivative of the partition func- 
tion with respect to 7 P . Such an antiferromagnet, which 
has no continuity with weak-coupling theory, is composed 
of composite excitations which can form excitonic bound 
states in the two-particle spectrum and hence is not in- 
consistent with the excitonic modes found in the mid- 
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infrared absorption of numerous parent cuprates 21 ^. In 
essence, the composite excitations described by the coef- 
ficients 7p and jp represent the orthogonal (they never 
lead to a turn-on of the spectral weight in the same en- 
ergy range) low-energy degrees modes that render the 
original UV problem weakly coupled. That such new 
degrees of freedom emerge as the dispersing modes is a 
typical feature of strong coupling. In fact, an analogy can 
be made here between our demonstration that the prop- 
agating modes at strong coupling in the Hubbard model 
are composite excitations (not electrons) mediated by an 
auxiliary field that has no bare dynamics with 't Hooft's 21 
demonstration that meson states, not free quarks, also 
mediated by a non-propagating auxiliary field, are the 



dispersing modes in QCD in 1 + 1 dimensions. Our anal- 
ysis suggests that a fixed point underlies the formation 
of such composite excitations. Whether the f3— function 
can be calculated within this formalism remains the out- 
standing question. 
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